Is Warm Inflation Possible? 
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I. INTRODUCTION 



First models of inflationary cosmology |l],|| were based on an assumption that inflation appears as a result of high- 
temperature phase transitions with supercooling in the early universe [||. This idea survived less than two years, after 
which it was replaced by chaotic inflation which does not require initial thermal equilibrium and phase transitions [|J . 
There were several reasons why this happened; here is a short list [H: 



• The assumption that the universe was hot from the very beginning is not necessary in the chaotic inflation 
scenario. Moreover, the assumption that inflation begins only after a long stage of cooling implies that the 
universe from the very beginning must be very large and homogeneous. This means that such models do not 

■ provide a complete solution to the homogeneity and flatness problems. 

on : 

• The theories where inflation occurs after supercooling require fine tuned effective potentials satisfying "thermal 
00 ■ constraints". 



The infiaton field typically has a very weak interaction with other fields, so it may be out of thermal equilibrium 
in the early universe. 



• Even if the infiaton field is in thermal equilibrium, it takes a lot of time for it to roll to the minimum of the 
temperature-dependent effective potential. In many cases, inflation occurs while the field rolls down, so when 

' it arrives to the minimum of the effective potential at <f> = 0, the temperature vanishes. 

• Even if the temperature does not completely vanish when the field falls to the minimum of the effective potential, 
it vanishes during supercooling, so the thermal effects become irrelevant for the description of the tunneling and 

$_i ' of the subsequent stage of slow rolling of the infiaton field 

For all of these reasons, the idea that inflation is somehow related to high-temperature effects was almost completely 
abandoned. However, recently there was an attempt to revive it in the context of the so-called warm inflation scenario 

§• 

The main idea of the scenario can be formulated as follows. Suppose the universe is in a state of thermal equilibrium, 
and the field <f> slowly rolls down to its minimum. When the universe expands, its temperature tends to decrease as 
a(i)" 1 , where a is the scale factor. Therefore one expects that the temperature in an inflationary universe falls down 
exponentially and immediately becomes irrelevant. However, if the scalar field interacts with other particles, it may 
transfer some of its energy to the thermal bath. This may keep the temperature from falling to zero. 

This regime may continue in a self-consistent way if the amount of particles produced due to the interaction of the 
scalar field (j> with thermal bath is large enough to keep the field from a rapid fall to the minimum of the effective 
potential. 

However, we immediately see a lot of problems here. During each time H^ 1 the universe expands e times, and 
the energy density of previously existed ultrarelativistic particles becomes e times smaller. Therefore most of the 
particles during warm inflation at any given moment should have been created during the previous time interval 
~ 0.2if _1 . If the energy release is going to keep the field from rapidly falling down, then the energy released each 
time 0.2H~ l must be very small. Otherwise the field rapidly falls down, and there is no inflation. But if the energy 
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release is small, then the total number of particles in the warm universe must be very small, and therefore their 
interaction with the scalar field <fi may be too small to keep it from rapid falling down. 

Despite this problem, the basic idea is rather interesting and it should not be discarded without a serious investi- 
gation. It would be very interesting to see, in particular, whether this idea could help to realize inflationary universe 
scenario in the context of string theory. 

The results of investigation of warm inflation indicated, as expected, that this regime is very hard to obtain ||. 
However, the methods used in this investigation were rather complex and not very intuitive. Therefore the basic 
mechanism of warm inflation remained obscure, and there remained a hope that one can obtain a good realization of 
this scenario by considering slightly more complicated theories of elementary particles. 

In this paper we will try to give a simple and intuitive explanation of the mechanism of dissipation of the energy of 
the inflaton field, which will help us to understand the main problems of the warm inflation scenario. Then we will 
confirm our expectations using more rigorous methods of non-equilibrium quantum statistics. 



II. INTUITIVE ARGUMENT 



In order to explain the basic idea of warm inflation, as well as its shortcomings, in this section we will give a simple 
and intuitive discussion of this scenario in the theory of a massive inflaton field <f> interacting with a massless field x : 

l = \{d»4>? + \{d, x f - - . (i) 

If one neglects interaction between these two fields, equation of motion for the scalar field <f) has a familiar form 

4> + 3H<j>= -m 2 cj). (2) 

The term 3H(f> represents energy loss of the field <p due to expansion of the universe. This term, which is similar to 
the viscosity term in an equation of motion of a pendulum in a viscous liquid, slows down the motion of the field 
4>. As a result, the potential energy of the field changes very slowly, which under some circumstances may lead to 
inflation ||. 

One could expect that interaction of the fields <f> and x ma y lead to an additional energy loss of the field 4> due to 
creation of x-particles (reheating), which can be represented by adding another friction term to this equation Jo|,|lO|] : 

4> + T4> + = -m 2 (t). (3) 

If T S> H, the field <j> will roll down more slowly than in the absence of interaction. Therefore one could expect 
that inflation may continue for a longer time |)|, and may occur even in such theories where otherwise it would be 
impossible M]. 

Further development of the theory of reheating has shown that the situation is more complicated. Addition of the 
term Tcb effectively describes energy loss due to particle creation only at the stage of oscillations of the scalar field <j> 
jiJI , and only in the case when these oscillations do not lead to parametric resonance [|l3| . It would be inappropriate 
to use this equation for the description of energy loss of the field <j) during inflation, assuming, as usual, that density 
of all particles during inflation is exponentially small. 

However, if the universe is hot, and there are many x-particles in the thermal bath, then in a certain approximation 
the motion of the scalar field (j> is indeed described by Eq. (^) . This effect is the basic ingredient of the warm inflation 
scenario. The standard description of this effect is very complicated because it involves summation of series of higher 
order diagrams in the non-equilibrium quantum statistics p^ , p^| . We will describe this method later. In this section 
we would like to give a simple interpretation of this effect which will simultaneously tell us whether this effect can be 
significant enough to give rise to an inflationary regime. 

Let us write equation of motion for the field <f> taking into account interaction between the scalar fields 4> and x va 
the Hartree approximation: 

4> + 3H4> + m 2 <f> + g 2 cl)(x 2 ) = 0. (4) 
In a state of thermal equilibrium one has |^,|^] 

(2tt) 3 J u x 
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where 

n x = — > = Vp 2 +9 2 ^ 2 ■ (6) 

A exp -if — 1 

We will assume that gcj) <C T, because in the opposite limit all thermal effects are exponentially small. In this case 

5 (>-!?)+- < 7 » 

Now suppose the scalar field 4> changes by A<f>. This changes masses of x-particles g<j), and the equilibrium value 
of n x q should change correspondingly. However, this change cannot happen instantaneously; it occurs with a time 
delay At ~ r~ , where T x is a decay width of x-particles at finite temperature. This means that n x deviates from 
its equilibrium value n e x q : when the scalar field reaches its value <f>, the occupation number remains equal to n x q at 

earlier time, when the field was equal to <j> - A<j> ~ <j> - <j)At, with At ~ T^ 1 0: 



dn eq 



-^r- 1 . (8) 



Therefore (x 2 ) in equation (^) slightly differs from {x 2 )eq- 

, 1 f n x (p)d3p [d 3 pdn x * 

Note that T" 1 depends on momentum p, it is small for p ^T. The leading contribution to the last integral is given 



by particles with typical momentum p ~ T. In this case, as we will show in Sect. HID 

g 4 T 



Evaluation of the last integral yields: 



r x ~ . (10) 



(x ) = (x U + ^ = n [i - - F ) + • (ii) 



where (7 = 0(10). 

With an account taken of this correction, Eq. (Eh can be represented in the following form jl6| : 

C6 2 C6 2 a 2 T 2 

<f>+^<f> + 3H0 + m 2 + g 2 <f>( X 2 )e q « + + 3fty + m 2 + ~^4> = 0. (12) 



In Sect. HID we will obtain this equation by a more rigorous, but somewhat more complicated method. The simple 



and intuitive approach which we have used in this section allows us to easily understand the nature of the viscosity 

term -f-<j) and evaluate its possible significance. 

Indeed, at the first glance this term could be very large. In fact, it is not suppressed by any powers of the small 
coupling constant g, and in certain cases -^f-<t> may be greater than the usual viscosity term 3H^>. This was the main 
reason why it was expected that this additional viscosity term could play an important role in inflationary theory. 

However, let us remember the procedure of the derivation of the term The leading term in the expression for 

(x 2 )eq m Eq- (@) is ^-independent: (x 2 )eq ~ T 2 /12. The term jjrfi> appears because the subleading term in Eq. (Q) 
depends on 4>, and because this term slightly differs from its thermally equilibrium value Thus the term ^-</> 

appears in Eq. ( |l2| ) as a correction to a correction, because <C <C j^. In other words, ^-(j) -C <? 2 0(x 2 )eg> 

so the new viscosity term in Eq. (|l^) is always much smaller than the usual thermal correction to the equation of 
motion of the field 0. This is an important fact, which will help us to evaluate plausibility of warm inflation. 

As an example, let us consider two limiting cases: 

Case I. gT < m. 

In this case m 2 3> g 2 (f>(x 2 )eq S> —$-4>, so the new viscosity term is completely irrelevant for the description of 
the evolution of the field </>, as is clear from the above intuitive derivation of the equation of motion. 
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Case II. gT > m. 

If g4> > T, then all thermal effects disappear. If g<f> < T, then rn 2 </> 2 < g 2 <f> 2 (\ 2 ) eg ~ g 2 (/) 2 T 2 < T 4 . In this 
case equation of state is determined by ultrarelativistic matter, p « p/3, and inflation cannot occur. 

In conclusion, new viscosity terms in the equation of motion of the inflaton field appear as a correction to a correction 
to the usual high temperature terms such as g 2 4>(x 2 }- These terms appear if one wants to write equation of motion 
of the scalar field in terms of the thermal equilibrium quantities such as (x 2 )eq, which differs only very slightly from 
its out-of-equilibrium value (x 2 )- If the leading thermal effects cannot result in existence of an inflationary regime of 
a new type, then the sub-subleading terms cannot do it as well. 

In this investigation we did not really evaluate the viscosity term —tf-<$>\ we only used our knowledge that this term 

appeared as a correction to a correction to the term j^T 2 , so it should be always much smaller than J2T 2 . It is very 
instructive to study the same issue by solving equation ( |l2|) directly, without making any a priori assumptions about 
the magnitude of the viscosity term. 

Suppose one has a stage of warm inflation in the regime where T > g<j), but m 2 3> g 2 T 2 (since otherwise one either 
does not have thermal corrections, or one has a noninflationary radiation dominated universe). In this regime, and 

assuming that ^- > H, Eq. (|12|) acquires the following form: 

0+^^ + 77^ = 0. (13) 

In the regime of slow rolling one can neglect the first term and obtain a solution 

,o 2m 2 Tt 



C 



(14) 



assuming, for simplicity, that the temperature remains constant. One can easily verify that for gm < T < and 



(j) 2 > Ejlf- this solution is compatible with the inequalities T > g<fi, but m 2 ^> g 2 T 2 , and the term <fi could indeed be 
neglected in Eq. (|l^) . Since it is not a usual oscillatory regime, one could think that we have an interesting candidate 
for the warm inflation regime. 

However, Eq. ( |l4|) tells us that the value of the field completely changes (vanishes) within the time 5t — 2m°T • 
As one can easily verify, in the regime which we study now (gift <C T and m 2 3> g 2 T 2 ) this time is much smaller than 
the typical decay time of the x particles: 

rst ~ « — « 1 . (15) 

x 3847rm 2 384vr V ' 

Thus the simple linear approximation (j) — A(f> ~ (f> — </>At, with At ~ r~ , which we used for the derivation of Eq. 
( |l2| ) does not work in application to our problem: The typical duration of the linear change of the field (f> is much 
shorter than T" 1 . This means that Eq. ( |l2| ) is incorrect in the only regime where warm inflation could appear in 
this model. In our previous qualitative analysis of warm inflation in this model we avoided this problem because 
we did not used the expression —f~4> for the viscosity term, which we obtained using the incorrect assumption that 

4> — A(f> ~ cf) — (/T^ 1 . Instead of that, we used the fact that the viscosity term is always much smaller than ^T 2 <p. 

To study motion of the field <fi in this model one should remember that the typical time when a significant change 
of the scalar field occurs (such as the oscillation time m -1 ) is much shorter than T x . Therefore the number of x 
particles practically does not change during the oscillations. To study the motion of the field <f> m this case one should 
evaluate (x 2 ) in the regime where n x remains constant, see e.g. |0. The behavior of the field <fi in this regime is 
completely different from what one could expect on the basis of Eq. ( |l^ ) . 

This example shows one of the many obstacles which make it so hard or perhaps even impossible to implement the 
idea of warm inflation. In order to verify our conclusions, we performed a more detailed investigation of this issue 
and studied several different models where warm inflation could occur. Before starting its detailed description we 
remind the reader that the above intuitive derivation of equation (|l^) is based on the following two most fundamental 
assumptions or conditions: 

(A) (f> interacts with particles whose effective mass is much smaller than the temperature. 

(B) <fr should not evolve significantly during the relaxation time, At, of such particles, namely, \<j)\ At <C <fi- (Otherwise 

the quasi-stationary inflationary expansion could not be sustained.) 
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A complete discussion of warm inflation is rather involved; in addition to these two conditions many other constraints 
have to be satisfied. However, as we will see in order to show that the scenario is not feasible in all models that we will 
consider it is enough to use only the simplest constraints based on the above two conditions. In the next section we 
review derivation of the equation of motion of the scalar field in more generic cases in terms of quantum field theory 
at finite temperature, where these two conditions (A) and (B) are also indispensable. It has been advocated in the 
literature J^] that this method is quite appropriate for application to warm inflation where the scalar field presumably 
moves slowly. 

As we have demonstrated above for a particular model of two scalar fields with interaction ^-</> 2 x 2 , however, 
does evolve significantly during the relaxation time even in its slow roll-over regime, so that the condition (B) is not 
satisfied. Put it differently, if we attempt to realize warm inflation consistently with the derivation of the equation of 
motion in the form of Eq. (|l^) , we inevitably find that the e-folding number of inflation is constrained to be much 
smaller than unity. In fact, a similar result has also been obtained in || for this particular model. However, there 
was a hope that this result is not generic, and several possible ways out to increase the number of e-folds have been 
proposed |^|]. In Sec. [V we demonstrate for several other models that the number of e-folds of inflation must be much 
smaller than unity. As will be seen there, this conclusion remains unchanged even when we take the proposals of Q 
into account. 



III. EFFECTIVE EQUATION OF MOTION IN THE THERMAL BATH 



As a preparation to quantitative discussion which confirms the above intuitive argument here we review a field 
theoretic approach to derive an effective equation of motion in the thermal bath. The standard quantum field theory, 
which is appropriate for evaluating the transition amplitude from an 'in' state to an 'out' state for some field operator 
O, (out \ O \ in), is not suitable to trace time evolution of an expectation value in a non-equilibrium system interacting 
with a thermal bath. For this purpose we should use the in-in formalism, which was introduced by Schwinger and 
developed by Keldysh jisfl . Here, extending Gleiser and Ramos []l5| and Yamaguchi and Yokoyama |^0[ , who followed 
Morikawa [|14| , we derive an effective Langevin-like equation for a coarse-grained field in the case the scalar field is 
interacting with not only other scalars but also fermions. 



A. Non-equilibrium quantum field theory 



Let us consider the following Lagrangian density of a singlet scalar field p interacting with another scalar field x 
and a fermion tjj for illustration. 



1 \2 1 „,2 2 



1 



-\<p 



1 



,2,, 2 



1 



X - 7.9 X V + #7^V - m^ipip - fipipip . 



(16) 



In order to follow the time development of ip, only the initial condition is fixed, and so the time contour in a 
generating functional starting from the infinite past must run to the infinite future without fixing the final condition 
and come back to the infinite past again. The generating functional is thus given by 



= Tr 



T ^exp 
T+ (exp 



i / ( Jp + Kx + r]ip + rjtp 



i / {J+ t P+ + K+x+ + V+i>+ + 



x T_ exp 



(J-P- + K-X- + V-i>- + V-^*-) 



(17) 



where the suffix c represents the closed time contour of integration and X + a field component X on the plus-branch 
(-co to +00), X_ that on the minus-branch (+00 to —00). The symbol T represents the time ordering according to the 
closed time contour, T + the ordinary time ordering, and T_ the anti-time ordering. J, K, 77, and fj imply the external 
fields for the scalar and the Dirac fields, respectively. In fact, each external field J + (K + , fj + ) and J_(i^_, rj-,fj_) 
is identical, but for technical reason we treat them different and set J + = J-(K + = K-,t] + = 7)_,?7 + — V-) on b7 a ^ 
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the end of calculation, p is the initial density matrix. Strictly speaking, we should couple the time development of the 
expectation value of the field with that of the density matrix, which is practically impossible. Accordingly we assume 
that deviation from the equilibrium is small and use the density matrix corresponding to the finite-temperature state. 
Then the generating functional is described by the path integral as 

Z[J,K, V ,f)}=exp(^iW[J,K,r ] ,fj}^ = Jv l pJv X Jvipjvre X piS[tp, X ,i,,j,J,K,r ] ,fj}, (18) 

where the classical action S is given by 

S[<p, X , 4>, $,J,K, V ,fl}= J d 4 x {£ + J{x)y{x) + K(x) x {x) + V (x)^(x) + fj(x)$(x)} . (19) 

As with the Euclidean-time formulation, the scalar field is still periodic and the Dirac field anti-periodic along the 
imaginary direction, now with <p(t, x) = ip{t — x), X (t, x) = x(t — i(3, x), and ip(t, x) = —ip{t — i(3, x). 
The effective action for the scalar field is defined by the connected generating functional as 

T[4>] =W[J,K,T],fj]- I d 4 xJ(x)<j)(x) , (20) 

J c 

where <j>(x) = SW[J, K,r],fj]/SJ(x). 

We give the finite temperature propagator before the perturbative expansion. For the closed path, the scalar 
propagator has four components. 

= / Tr[T +X (x) X (x')p] Tr[ X (x') X (x)p] \ ( 
~ { Tr[ X (x)x(x')p] Tr[T_ X (x) X (x')p]J [ZZ) 

Similar formulae apply for ip field as well. 
Also, for a Dirac fermion we find 

Tr[T + iP(x)$(x')p] Tr[-i>{x')i>{x)p] \ , , 

Tr[ip(x)ip(x')p] Tr[T_ip{x)il>{x')p)) 1 ' 



B. One-loop finite temperature effective action 

The perturbative loop expansion for the effective action T can be obtained by transforming p — > ipa + £ where <po 
is the field configuration which extremizes the classical action S[ip,J] and £ is small perturbation around po. Up to 
one loop order and 0(A 2 , g 4 , f 2 ), depicted in Fig. effective action T becomes 

r[0 c ,<M] = J d i x^x)[-n-M 2 ]Mx)-^(HA(x)4(x) + ct )c (x)ct ) 3 A (x)) 

-J d 4 x J d 4 xAi(x - x') 
-2 J d A x J d 4 x'A 2 (x - x')4> A {x)<l> c {x') 



&{x)(t> c {x)(t> 2 c {x') + -4>a{x)4> c {x)4>\{x') 



+ \jd A x J d 4 x' 



Bi(x - x')(j) A (x)(j)A(x')(f)c(x)(j) c (x') + B 2 (x - x')4>a(x)4>a(x') 



(25) 



where 
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M 1 =m 2 + g' 



2 f d 3 q l + 2n x (q) X f d 3 q 1 + 2n v (q) 



(2tt) 3 2u: x {q) 2 J (2tt) 3 2u v {q) ' 



Ai (x - x') = 2 ff 4 Im [G*(x- x') 2 ] 6{t - t') + —Im [G*{x- x'f ] 6{t - t') 
A 2 (x - x 1 ) = / 2 Im [ S^(x - x')S p F a (x' - x) ] 0(t - i 7 ) . 

Bi(x — x) = 2g 4 Re [G x (x — x') 2 ] + -yRe [G^(x — x') 2 ] . 
B 2 (x-x') = -f 2 Re[S^ p (x-x')S^ a {x'-x)] . 



(26) 
(27) 

(28) 

(29) 
(30) 

(31) 
(32) 



The last term of (|2^) gives the imaginary contribution to the effective action T. We can attribute these imaginary 
terms to the functional integrals over real auxiliary fields £i(x) and £2 (x) H] to rewrite (Eq) as 



exp(*r[0 c ,0 A ]) = Jvfr Jv^Pi[^}P^ 2 ] exp{^ e// [0 c ,0 A ,6,6]} 



Se&[(l>c,(l>A,£i,&] = ReF + / d 4 a;[^i(a;)^ c (a;)0A (») + £,2(x)(/)a(x)} 



where 



Here £1 (x) and £2 (x) are random Gaussian fields with the probability distribution functional 



Pilti] =Mexp 



~ / d 4 x / dV^Wi?- 1 ^-^)^^') 



(i = 1,2). 



(33) 



(34) 



(35) 



where A/i is a normalization factor. They induce random noise terms in the effective equation of motion of <f> as a 
result of the interactions with the thermal bath. 



C. Equation of motion 

Applying the variational principle to S e g, we obtain the equation of motion for 

(Steffi 0c, <M, 



. 



(36) 



From fl3j) and (E5J), it reads 



(□ + M 2 )0 c (*) + -^(x) 



and 



+ c (a;) / dV / dt'Ai(z-<)<^(<) + 2 / dV / di%(:r - x')cj> c {x') 

) J J — OG 

= <£c(a)£i(a;) +6 (x) , 
(ti(x)titf)) =B i {x-x'). 



(37) 



(38) 



Though A\ and i?i has two contributions from \ an d V? fields, they have the same properties except for the values of 
coefficients and masses. For the moment, we consider only the contribution from (p field for simplicity and omit the 
suffix c. The right hand side of (p7|) are the noise terms, while the last two terms of the left hand side are combination 
of dissipation term and one-loop correction to the classical equation of motion which would reduce to a part of the 
derivative of the effective potential, V^ s (4>), if we restricted <j>(x') to be a constant in space and time. 

The above equation ( |37| ) is an extension of equation (3.2) of Gleiser and Ramos Q in that we have incorporated not 
only self-interaction but also interactions with a boson \ and a fermion ip. ft is nonlocal in space and time. The spatial 
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nonlocality does not bring any difficulty here since the scalar field is presumably nearly homogeneous during inflation. 
So we only need to consider contributions with zero external momentum, that is, we can put <f) c (x',t') — <p c (x,t') in 
the integrand. 

With this approximation the correlation function of the bosonic noise ( (3l|) with g — 0, for example, becomes 

= ^S a (x-x') J -^Re[G$(g,t-f)] 2 . (39) 

We thus obtain spatially uncorrelated noise. For the scalar field averaged over a volume V, the amplitude of noise 
reduces in proportion to V~ x l 2 , which implies that noise term should be omitted in the equation of motion of spatially- 
averaged homogeneous field. On the other hand, the temporal nonlocality is very important in deriving the dissipation 
term as seen below. 



D. Dissipation term 



The equation of motion ( p7\ ) derived above has contributions representing the dissipative effect in the last two terms 
of the left hand side. Since these terms are nonlocal in time, in order to extract local terms proportional to <j> one 
should assume that the field changes adiabatically |uj |L5j , or put 

4> n (x', t') ~ (/> n (x', t) + n(t' - t)<f> n -\x', t)4>(x', t) (40) 
in the integrand of (^). Then these terms will read 

4> I d 3 x' [ dt'A 1 {x-x')(/) 2 (x') + 2 [ d 3 x' I dt'A 2 (x-x')4>(x') 

J J — oo J J — oo 

= 4> 3 {t) I d 3 x' I dt'A 1 {x-x') + 2(f>(t) I d 3 x' I dt'A 2 (x-x') 

J J — oo J J —oo 

+20 2 (i)0(i) / d 3 x' [ dt'(t' -t)Ax(x-x')+2<j>(t) [ d 3 x' [ dt'(t' - t)A 2 (x - x'). (41) 



The last two terms are dissipation terms. They would vanish if we used bare propagators, as a manifestation of the 
fact that the dissipative effect is intrinsically a non-perturbative phenomenon and cannot be investigated from the 
perturbation theory (2lJ . In order to obtain a finite result we should use full "dressed" propagators instead |l4|,[l5|,|| . 

The viscosity from scalar interactions have been fully investigated in [|l5| and || , so we simply quote their results 
here. The full propagator of cj> reads 



= — /[l + n B (oj$ - zr )] e -*W-«V)|t-*'| + ns(cj0 + ir ^ e <(**+<r,)|i-t'| | j (42) 



(2tt) 3 

where = (e^"* — l) 1 and ui^ = yj k 2 + irL^ T with to^t being the finite-temperature effective mass of tp. Here 

is the decay width related with the imaginary part of the self energy, E, as 

Im£ ^ A 2 T 2 

in the limit T 3> m^x- We find the contribution of the self- interaction to the viscosity 

AV(W) / d 3 x> f dt>(t'-t)lm[G$(x-xr]^ ^(3 f ^^i±^ (44) 

96 , / T 
~ — m 



7rT \ m^T 
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in the high temperature limit and in the case g = / = 0. Note that the above expression was first obtained by Hosoya 
and Sakagami [jl6|| by a different method which is intuitively more appealing as discussed in Sec. Ejl. 

The contribution from the interaction with \ can be obtained similarly. Assuming \ has no self-interaction its 
width is given by 



r y ~ ^ (45) 

x 1927T 



in the high-temperature limit and the relevant part of the viscosity terms is given by 

4 5 V(*M*) J d 3 x' J dt'(t' -t)lm[GF(x-x') 2 ]~^ln(J^j<j> 2 j>. (46) 

In both cases the viscosity terms due to scalar interactions are of the form 4> 2 4>/T ||. 

On the other hand, the dissipation due to Yukawa interaction is calculated in the Appendix. In the high temperature 
limit T 3> m^T, rn^x we find 

2f4> f d 3 x' f dt'(t' - t)lm[S^(x - x')S^ a (x' - x) } ~ ^C(3)T0 ~ 

J J — oo 

In both cases the viscosity coefficients would be exponentially suppressed and would not play any important role if 
the high-temperature conditions were not satisfied. 

IV. FEASIBILITY OF WARM INFLATION 

In this section we study if the warm inflation driven by the viscosity term is possible in the case the viscosity term in 
the effective equation of motion is most effective, namely, in the high-temperature limit when the viscosity coefficient 
is given either by ( fiil ) or (ffTj ) depending on the interaction. We consider these cases separately for both new Q and 
chaotic [jlj inflation scenarios. 

Before analyzing specific models we formulate generic conditions to satisfy. We are interested in the new possibility 
that the slow-rollover inflation is realized due to the thermal viscosity (Q or (f47|), so the effective equation of motion 
should read 

C v ^ = -V'[(j>], (47) 

that is, we require C v 3> 3H. The inflaton's energy released through this viscosity term presumably goes to radiation, 
whose energy density, p r , satisfies 

^L = -4H Pr + Cj 2 . (48) 
at 

Since C' v strongly depends on the radiation temperature, p r should not change too rapidly in time in order to sustain 
quasi-stationary inflation. So the creation term in ( p8| ) should balance the redshift term. As a result we find 

Pr = ^T^^, (49) 
' 30 4H ' V ' 

which gives the radiation temperature as a function of <j). In contrast to |S, where the temperature has been fixed 



to its initial value, we perform a consistent analysis by using the value of T calculated from (49), which would help 
to increase the viscosity from bosonic interactions || . Here is the total effective number of relativistic degrees of 
freedom. We normalize it by 150 and denote g»/150 = g*N below. In order that the universe is inflating the potential 
energy density should dominate over p r and the kinetic energy density. 

Finally a number of conditions must be satisfied to justify the derivation of the effective equation of motion. In 
order that the radiation created from the inflaton thermalizes sufficiently rapidly, we need 

depending on the nature of interaction. In addition, the adiabatic conditions 







(51) 



must be fulfilled so that the viscosity term is proportional to <p. The inequalities (^0|) and ([31]) are essential to realize 
quasi-stationary state of inflationary expansion. Finally, but most importantly, the high-temperature conditions 
T 3> rn^T: m xTi m-ipT must also be satisfied, otherwise the viscosity term would be exponentially small and our 
entire discussion would break down. 

Since the failure of the Case II in Sec. [n] is evident, we here consider the Case I, where finite-temperature 
correction to the effective potential is sub-dominant but still the viscosity term appears sizable at first glance (until 
we convince ourselves it is not, through the intuitive argument in Sec. [H].) Thus we can make the following list of the 
inequalities to be satisfied to realize the desired scenario. 

(i) c v > m. 

(ii) v[4>] » \4> 2 . 

(iii) V[4>] > pr. 

(iv) r> > H, T x > H, or I> > H. 



(v) r » 



r x » 



, or > 



(vi) T > m 0T . 



(vii) T > m X T or T > m^r- 

(viii) Finite-temperature correction to the effective potential being sub-dominant. 

Inequalities (v) are nothing but the condition (B) in Sec. O, and (vi) and (vii) stand for the condition (A). As will be 
seen below, we can essentially rule out all the models we consider only in terms of these two conditions. 

In addition to the constraints listed above, there is a constraint following from the investigation of density pertur- 
bations produced during warm inflation. 

The standard expression for density perturbations produced during inflation in the cold-matter dominated universe 



Sp _ 6 BM 
P 5 4> 

If one uses the standard estimate for inflationary perturbations 



H_ 

2tv 



one gets 



5_p 
P 



6 H 2 
5 2tt0 



(52) 



(53) 



The first of these two equations holds for the warm inflation as well. However, the second one should be modified 
because the amplitude of fluctuations during warm inflation is greater than -M- for two different reasons. 

First of all, the wavelength of perturbations of a field with mass m 2 <C H 2 which freeze during inflation usually is 
0(-ff _1 ) because of the friction term 3£f0. This allows one to make the standard estimate 8(f) ~ S- by calculating 
the amplitude of vacuum fluctuations with the wavelength 0(-ff _1 ). However, during warm inflation the amplitude 

of perturbations of scalar fields with momenta k -C T is enhanced by the factor ~ J ? because of the additional 



contribution of thermal fluctuations. This leads to an estimate 8cj) ~ y P^ j- 

But this is not the only effect which should be taken into account. Indeed, during warm inflation the friction term 
is [C v + 3H)(f>, and it is assumed that C v ^> 3iJ. As a result, it may happen that the fluctuations of the scalar field 
freeze not when their wavelength approaches H^ 1 , but much earlier, when their amplitude can be much greater. This 
may lead to an additional increase of the magnitude of perturbations of scalar fields produced during inflation. 

We will not perform here a detailed evaluation of density perturbations in warm inflation because we did not find 
any model where this scenario can be realized. Indeed, we will be able to rule out warm inflation in all models 
which we will consider even without using the theory of density perturbations. However, one should keep in mind the 
necessity to study constraints based on the theory of density perturbations, because usually these constraints lead to 
the strongest restrictions on the structure of inflationary models. 
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A. Chaotic inflation with viscosity from bosonic interaction 



Here we consider chaotic- type inflation with a potential V[(f>] = -^(j) 4 and the viscosity arising from self- interaction or 
other bosonic interactions proportional to (f> 2 , such as ^g 2 (f> 2 x 2 ■ In the usual chaotic inflation slow-roll of the inflaton 
is realized due to the Hubble friction and it is effective only when <\> > O.SMpi, but if thermal viscosity discussed 
above is effective, we might have inflation with much smaller <j). 

In this case the viscosity coefficient and the Hubble parameter are given, respectively, by 



96c0 
ttT 1 



H 



VttA 
3Mpi ' 



(54) 



Here = ln(T ' /tti^t) if </> has self-interaction only. But can be increased if it interacts with additional scalar fields 
Xi's through the interaction 



El 



2 9j(P Xy 



(55) 



Then C v is given by 



96 

7rT 



A 2 



a 2 + E^ 



111 



T \ 



T 



m XjT 



(56) 



where rn Xj T <C T is the finite-temperature effective mass of Xj field H]. For simplicity, we identify C0(3> 1) with the 
number of \j fields, which corresponds to the case \n{T/m Xj T) = 2. We also assume all the coupling constants gj 
take the same value Qj = g. 

From ( f47j ) and (|4^) the temperature is given by 



T = 1.6x lO-'gJc^XiM*^* 



(57) 



We can then convert the inequalities (i)-(vi) to the constraints on the range of <f> and on other model parameters as 
follows. 

(i) ^«3.4x 10 4 <&v4A-*M PI . 

(ii) 4> » 2.9 x 10- u g-^X 3 Mpi. 



(iii) 0» 1.3 x lO-^c^Af M, 



pi- 



(iva) </> « 1.2 x lO- 4 gJc^\?Mpi. 

_1 __ 



(ivb) 0«5.5x W^gjc *g 3 . 



(va) A> 26c^ 1 . 
(vb) Xg- 4 < O.3Oc 

(vi) 0«1.2xlO- 5 .g^ 1 M Pi . 

(vii) « 4.3 x W-^g-^Xig- 3 . 

(viii) A > 3.8 x 1CT% 5 4 . 

The last condition is from the requirement that radiative correction due to \ does not change A. 
The e-folding number of warm inflation, if any, is calculated as 



N : 



<t>i 



H~T = 5-1 x lO^vf A-^ (<, 



(58) 
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where <f>i and cf>t are upper and lower bounds of <f> that satisfy all the above inequalities. 
Using (vii) and (vb) in ([38]) we find 

N « 3.6 x 10- V^.<T 4 « 1.1 x lQ~ A g-^. (59) 

Since x' s a l so contribute to g*, g~ N c<f, cannot be larger than 150 and it is apparent that N cannot even exceed 
unity. Hence we cannot realize warm inflation in this model. This conclusion is independent of the simplification 
\u{T/m X T) = 2, for the above constraint on N is primarily due to the condition (vii) or T> ti x t, and the fact that 
this condition is hardly satisfied implies we have overestimated the duration of warm inflation. 

Next for completeness we consider the potential V[4>] = ^m 2 ^ 2 + ji4> 4 with m 2 ^> X(f> 2 /12. In this case the 
temperature is given by 

T = 4.3 x lO-V*<V*mM*,0-*, (60) 

and the inequalities read 

(i) 4> » 2.8 x lO-Vjc^miMj,. 

(ii) 4> » 6.0 x 10- 2 g J c^mUl^. 

(iii) (j> > 9.1 x lQ- 2 g~j^c^miMl v 
(iva) « 9.5 x 10~ 5 g J c^XiMpi. 
(ivb) 4> « 4.5 x l^g-Jc^gHlpL 

(va) 0> 13c <f) 2 \- 1 m. 

_ i 

(vb) 4> > 4.4c^ 2 g~ 2 m. 

(vi) (f> « 8.0 x lO-^c^Mp;. 

(vii) « 9.4 x lQ- 2 g'Jc^g-imiM~ pi . 

(viii) < \/l2\~^m. 

The number of e-folds of inflation is formally given by 

N = J H-4- = 4.4 x 10 2 m- 2 Af p; 3 - <f> f * j . (61) 

From (vb) and (vii) we find m « 2.1 x IQ~ 7 g'^c^g 5 M Pi . Then inequality (vii) reads (f> <C 9.2 X lO^ 7 g~^c^g 3 M Pi . 
Using these inequalities in (|6l]) we find 

7V«7.3x IQ^g-Jc^. (62) 

Thus no matter how many scalar fields are interacting with the inflaton we cannot find warm inflation solution. Note 
that the above constraint on the e-folding number primarily comes from the condition (vii) orT> rn X T ■ This means 
that if we had used the correct value of hx{T/m X T) instead of putting it to be 2, the viscosity term would have been 
smaller and warm inflation would have been even more unlikely. So we can justify our simplification. 

Therefore the conclusion that the number of e-folds of warm inflation with chaotic-type potentials is constrained to 
be much smaller than unity || remains unchanged even when we use the consistent value of the cosmic temperature 
obtained from (|4^) , rather than fixing it to its initial value || . 

Note that these considerations look different from the arguments used in the Sec. |l[ but they lead to the same 
final conclusion. Now let us see how one can reach the same conclusions directly, using the arguments of Sec. |fi|. We 
will assume for simplicity that the field <fi interacts with one field %, and A <C g 2 ■ In this theory Eq. ([l2]) looks as 
follows: 

4> + C v 4> + 3H<j> + m 2 <f> + \<p 3 + g 2 <t>{x 2 ) eq - 0. (63) 
Let us consider two limiting cases: 
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Case I. g 2 T 2 < max {m 2 , Ac/> 2 }. 

In this case m 2 <j)+X(j) 3 3> g 2 (f>(x 2 )eq ^> C v <p, so the new viscosity term is completely irrelevant for the description 
of the evolution of the field (j), as is clear from the intuitive derivation of the equation of motion in Sec. E. 

Case II. g 2 T 2 > max {to 2 , X(j) 2 }. 

If gcf> > T, then all thermal effects disappear. If g<j) < T, then m 2 < ? !) 2 /2 + A<?(> 4 /4 « <?V(x 2 } eg ~ ffV? 12 <C T 4 . 
In this case equation of state is determined by ultrarelativistic matter, p w p/3, and inflation cannot occur. 



B. Chaotic inflation with viscosity from Yukawa interaction 



We next study the case the dominant contribution of viscosity arises from Yukawa interaction, so that we find 
C v ~ llc^T where = 1 in the case <f> has a Yukawa coupling to one species of Dirac fermion, but it can be larger 
if 4> interacts with more fermions. We identify c^, with the number of fermion species interacting with (f> with the 
universal coupling strength /. 

First we consider the case inflation is driven by a quartic potential V[<fi] — ji<fi 4 - The temperature is given by 

r = O.12 5 ; j *^*A*0*M* I . (64) 
The following inequalities must be simultaneously satisfied. 

(i) 0«O.78 3 JcfA-Uf Pi . 

(ii) 0«62 3 ;jc 4 A- 1 M Pi . 

(iii) 0>O.15 5 J Jv c^ 1 A3M P ;. 

(iv) 0«O.63 ff J A ,c|/ 5 A- 1 M Pi . 

(v) / 4 » 0.31c^A. 

(vi) « 1.4 x lQ-±g-^ l \- x M Pl . 

(vii) « 2.5 x W^g-^f^XiMpi. 

(viii) 0»1.6x lO-^clpX-K 

The last inequality comes from the condition that the finite-temperature correction to the effective potential is small, 
that is, c^/ 2 T 2 /6 < X(j) 2 /2. 

The number of e-folds of inflation is expressed by 

N = J ' H d -i = 5.8 5 J4 A "* M p* " 4) ■ ( 65 ) 

Using (v) in (vii) we find <fi ^ 1-1 x 10~ 4 <7~nC^ A* Afp;. I nser t m g this limit to <f>i in ( |65j ) we obtain an upper bound on 
N as N -C 3.9 x 10 _3 <7~^-c,/,. Since g~^c^p cannot be larger than 300/7 we can conclude warm inflation is impossible 
here. In fact, we can explicitly show that there is no open parameter region that satisfy all the inequalities even in 
the case g~^c^, is maximal. 

Next we analyze inflation driven by the mass term V[<j>] = ^m 2 (j> 2 , in which the temperature is given by 

T = 0.l9g;*c-*mi M|^i (66) 
The following inequalities must be satisfied for successful warm inflation. 

(i) 4> <0.26 5 ;J c^m-iM^. 

(ii) <p > 0.025 5 * A rc^ 4 TO 2 Af p; 1 . 
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iii) 0> 0.170V* *to*M»,. 

iv) 0»13 5 ;jc7/- 5 m 2 M Pi 1 . 

v) > 1.4c^/ _2 m. 

vi) > 4.2 x 10 3 g* N c^m 2 Mpl. 

vii) « o.i3 5 ;j c ; 4 /-*miM| r 

viii) < 3.6 x 10 5 5*jvc^/~ 5 w 2 Mpj 1 . 

The last inequality is from the requirement m 2 3> ^-f 2 T 2 . 

It is not impossible to find the values of model parameters that satisfy all the above inequalities. This does not 
imply, however, that warm inflation is feasible in this model because the e-folding number of inflation, 



r<t>t i 

Y : I H- 

<j>i 



3.6 ff Jcjm-8M P * (cj>f - 4>j ) , (67) 



turns out to be smaller than unity as seen below, where 4>i and 4>f are upper and lower bounds satisfying all the 
inequalities (i) through (viii) as before. 

Inserting the upper bound (vii) to fa in (pij), we find 



N « Q3lgl N clr*m*M P ?. (68) 



From the consistency between (v) and (vii) we find m -C 7.4 x 10~ 5 g~^Cip f 3 Mpi, which, together with (68), implies 
N <C 2.7 x 10~ 3 g~ ] t I c t /j. Since g~^c^ cannot exceed 300/7 we can conclude warm inflation is impossible in this model, 
too. 

C. New inflation 

Next we consider new inflation driven by the potential 

fm= K 02 ~x) 2 - (69) 

Inflation is possible only for 

4><^\-^m. (70) 
First we analyze the case viscosity is dominated by bosonic interaction. The temperature is constant in this case: 

T = 4.8 x 1(T 2 g^c^X^mHl^. (71) 

We find the following inequalities. 

1 2 



(i) ^.»8.3x 10- 2 g J c^X-emzMpf. 



(ii) <p > 2.2 x 10- 3 ff ^c^AlmiAf 



(iii) m > 6.0 x IQ-ZgJc^XiMpi 



1 1 

PI- 



(iva) m < 1.3 x 10~ 4 g^ c 4 A 2 M 

i _ i 



(ivb) m « 6.4 x l(r 4 c, *gJ\*g 6 M Pl 



4 \± 3 
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(va) 0> 13c. 3 A" 1 ™. 

_i 

(vb) 4> > 4.4c^ 2 g- 2 m. 

(vi) m« 1.1 x 10-V^c^AiMK. 

—1 —1 1 2 - 

(vii) </> < 4.8 x 10~ 2 g^c^ 3 A« g~ x m^ Mp v 
(viiia) to > 1.4 x lO^^cT^Mpj. 

(viiib) to » 2.7 x 10- 6 9 ^4f 3M «- 

Here the conditions (viii) are required to ensure the symmetry remains broken and warrant the use of the zero- 
temperature potential, that is, (viiia) comes from the condition to 2 ^> AT 2 /4 and (viiib) from to 2 ^> c^g 2 T 2 /12. 
The number of e-folds of inflation is expressed as 

TV = 4.6 x lO^cf A ^ m ~ §M p/ (<f> 2 - 2 ) , (72) 

where <f>i and 4>j are now the lower and upper bounds of 4> which satisfy all the above inequalities, respectively. 

Inserting (vii) to (|7|) we find TV <C l.lg^^ X^ g~ 2 mi M pl 3 . The consistency between (vb) and (vii) sets an upper 

bound on m as m « 1.3 x 10~ 6 g~^c^ 1 X^g 3 Mn. These two inequalities implies N <C 1.3 x 10 -4 g~-A^; thus warm 
inflation is not feasible. 

So far we have used only inequalities (v) and (vii) as promised, apart from the generic condition (|7^) for new 
inflation. If we use other inequalities in addition, we can completely close the allowed region of the parameter space 
as follows. The consistency between ( |70| ) and (ii) sets a lower bound on to as 

to » 1.1 x 10- s g^clxHd Pl . (73) 

From (vi) and ( |73[ ) we find A < 22c! 1 . On the other hand, from @ and (v) wc find A > 1.7 x 10 2 cl x . Thus there 
is no allowed region for A to realize warm inflation consistently. 

Next we move on to the case viscosity is dominated by Yukawa interaction, when the temperature is given by 

T = 0.20g^c~ i \^mi(t>iMl, l . (74) 
The following inequalities must be satisfied. 

(i) 0»5.5.gJ Ar c^ 2 A-lTO 4 M P f. 

(ii) 2.1^c|A-imiMj r 

(iii) 0<2.1 ff Jc|A-ImlMp*. 

_ 1—1 1 5 3 —I 

(iv) <j> > S.Og^c^ 8 X~- f~2m2 M pl 2 . 

(v) 0> lAc^f~ 2 m. 

i i _1 3 -i 

(vi) cf> 56g 2 N c^X ^m^M Pl 2 . 

(vii) 4> < 6.8 x 10~ 2 g~*c~ V~^m§M| r 
(viiia) <j> < 3.2 x IO^cJa- ir^M"*. 
(viiib) </> < 5.3 x lO^VTV-lA-SmlMp;'. 
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(viiia) comes from the condition m 2 3> XT 2 /A and (viiib) from m 2 3> c^f 2 T 2 /6. 

In this case, contrary to the case the viscosity arises from bosonic interactions, there exists some allowed region in 
the parameter space, but the number of e-folds of inflation, 

N = jf ' jd0 = S.fjgjcf A-i m^Mpf (fa) - 0,f) , (75) 

turns out to be much smaller than unity as seen below. Here fa and fat are lower and upper bounds on <p obtained 
from the inequalities (i) through (viiib) and (|70| ) as before. 
From (vii) we find 

N « 2.1g-J c|A-i/-i m§ M~f . (76) 

Consistency between (v) and (vii) imposes an upper bound on m as m « 1.2 x 10~ i g~ ] ^c^A^ fMpi. Inserting it to 
© we find N < 6.6 x 10 3 g^^c^. Thus N is much smaller than unity no matter how many fermions are interacting 
with the inflaton. 

D. Shifted field model 

So far we have studied the possibilities of chaotic inflation and new inflation models driven by thermal viscosity 
term. As a result we have shown that none of the above models can accommodate inflationary expansion which lasts 
more than one e-fold. In fact, we have been able to rule out them simply from the high-temperature condition (vii) 
and the adiabatic condition (v). To cure this problem another model has been suggested in Q in which the scalar 
field x has the interaction term 



£ int = -E^ 2 (^-M) 2 X 2 . (77) 

j 

We are not aware of any particle physics motivation for this kind of interaction. Nevertheless it may be worthwhile 
to study this model because it helps to reduce the effective mass of \ when <j> is large and close to M. Thus one could 
hope that one may relax a constraint from the high-temperature condition (vii) and find a sensible solution for warm 
inflation g. 

The situation, however, is not that simple, as one immediately recognizes once he writes down the effective equation 
of motion in this model. Indeed instead of eq. (jl|), with C — Q&c^/tt, we find 

4> + ^r(0 - M) 2 cf> + <3H4> + m 2 + g 2 ^ - M) J^ix ))e q - 0, (78) 

3 

Thus the viscosity term vanishes at <fi = M, which makes the dissipation inefficient in the region <j) « M, 

Nevertheless one may still want to consider the possibility of warm inflation in the vicinity of <j) = M, where the 

high temperature condition, m x = g\fa— M\ <C T, is satisfied but the slight deviation from = M makes the viscosity 

nonvanishing. Below we consider this possibility in chaotic inflation driven by the mass term for illustration and work 

out inequalities to be satisfied as we did above. 

In this model the effective potential of (j> is given by 

V eS [fa = \m 2 fa + ^fr 2 (0 - M) 2 + ... (79) 
in the high-temperature limit. We must treat the cases with c^g 2 T 2 -C m 2 and with c^g 2 T 2 ^> m 2 separately. 

1. The case c^T 2 < m 2 . 

First we consider the case c^g 2 T 2 <C m 2 . Since we are interested in the regime <f> > \4> — M\, the second term in 
( [79| ) is entirely negligible in the inflaton's dynamics. The effective equation of motion in the slow-roll regime is given 

by 



1 
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CA = ^{<l>-Mfj> = -m% (80) 



and the radiation temperature is calculated from (|49|) as 

T = 4.3x lO^c^g'Jm^^-MyiM^. (81) 

Apparently it is divergent at = M . However, since = M is out of the slow-roll regime and will be automatically 
excluded using the inequalities below, this does not cause any problem. 
Now we list the required inequalities relevant to \\ 

(i) |0 - M| » Q.nc^gJmUHlp?. 



(ii) |0-M| >8.5 x lir'i-c^gJmt^M* 



PV 

i „ , i 

-2, 



(iii) |0- M| > 5.0 x 10" g t ^m"^4JW^. 

(iv) |0- M| « 2.1 x lO- 7 c7. g ;j.g 6 0- 1 Af2 ; . 

(v) |0 -M| >2.5 5 - 2 c^m. 

(vii) |0 - M| « 0.15c7<? Jm§0"Uf) r 
The number of e-folds of inflation is given by 



N 



l0 ~ M|/ d\d>-M\ -ii,-*!/ ii ii \ 

tf^U ^4xl0 2 c| 5 >- 2 M P! 3 f3 |0 -M|/ -|0 -M]/), (82) 

where the integration has been done near |0 — M\ assuming ~ M. Inserting (vii) to |0 — M|$ in ( p2[) one finds 

AT<0.4c~ f 9 J™MM;*. (83) 
From the consistency between (v) and (vii) it follows that 

m « 8.8 x lO-*cJ^ff B 0^M|,. (84) 

Using it in (]83) one finds that 

N « 9 x 10- 2 ^ 5 ;j#Uf P *. (85) 

This means that unless ~ M 3> Afp; one cannot obtain sufficiently long period of warm inflation. But then it 
is not interesting because for 3> Mpi we can realize chaotic inflation without the help of thermal viscosity. This 
conclusion follows from (v) and (vii). 

But if we consider other inequalities the situation becomes even worse. From the consistency between (i) and (iv) 
we find 

m0 2 « 1.3 x 10- 8 c-^gJg 8 M 3 Pl . (86) 

Inserting it to ( |83|) we obtain 

N^OMc^gJgi. (87) 
Thus warm inflation is ruled out for perturbatively meaningful values of the coupling constant. 
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2. The case c^,g 2 T 2 > m 2 . 



Next we consider the opposite limit, c^,g 2 T 2 3> m 2 . If g\(f> — M\ ^> T, then all thermal effects disappear. If 
g\4> — M\ <C T, one may still neglect contribution of the second term in ([79]) to the total energy density of the 
universe, because otherwise the universe would be dominated by radiation and there will be no inflation. This does 
not mean, however, that this term does not affect the motion of 0. On the contrary, the potential force is dominated 
by its derivative. As a result the minimum of the potential is shifted from the origin to the vicinity of = M: 

M 12m 2 M , . 

m = — 12m3 =M--_ -. (88) 

One may encounter two possible regimes by comparing the time scale of the friction, Tf = C" 1 oc T(0 — M)~ 2 , 
with that of 0's oscillation, r c = y/l2/ (c%gT). If the deviation from = M is sufficiently large to warrant ry < t d , 
the field is in the slow roll-over regime toward the minimum ( |8^ ) and warm inflation might be possible If Ty > t d , 
on the other hand, sits in the minimum ( |88| ) which is time-dependent through the temperature. Let us consider 
these two regimes in turn. 

First, if is rolling toward the minimum, we find the following effective equation of motion, 

^(0 - Af) 2 = -m 2 - ^.9 2 T 2 (0 - M). (89) 

The second term in the right-hand-side (RHS) dominates the potential force except in the close vicinity of = M, 
and we neglect the first term in RHS. In the r egime where the first term dominates over the second term we would 



come to the same conclusion as in Sec. IV D 1 



From ([49j) one can find the temperature 

T = 1.8 X lO^g^NC^m^Mpl. (90) 
We write down the required inequalities relevant to \- 



(i) |0 - M\ » 6.0 x lO^.V^m^M 



-l 
pi ■ 



(ii) \<f>-M\ » 1.6 x W & gl N c-^g 2 m 2 <t> 2 Mpl 

(iii) « 3.1 x Ur^g-Jclm-Hlpp 

(iv) .g 4 » 3.3 x 10-% 5 ^. 

(v) |0 - M| » 5.3 x lO^V^m 2 ^ 2 . 
(vii) |0 - M\ < 1.8 x lO^g^NC^g^mcjjMpl. 

The number of e-folds of inflation is calculated as 

N ee P*~ M[ ' H^^l = 6.4 x l0-^gJ4m- 2 ^ 2 M Pl (|0 - M\ 2 - |0 - M\ 2 ) , (91) 

J\<i>-M\i 



where the integration has been done over |0— M\ assuming ~ M, Inserting (vii) to (91 ) we find N <C 2.1 x 10~ 4 g~^c^. 
Hence this regime does not lead to inflation. 

Next we suppose that has fallen to the minimum (^8|) . Since the temperature is presumably gradually decreasing 
during the warm inflation, changes with time according to 

• 24m 2 M ■ , . 

-T. (92) 



c<pg 



2^3 



1 Since we are interested in the feasibility of viscosity-driven warm inflation we assume C v > 3H and hence the only relevant 
time scale of friction is Tf = CjT . 
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We would like to see whether this time- variation of <fi may create a sufficient amount of radiation by the energy release 
through the viscosity term to support quasi-stationary stage of warm inflation. 
If such a stage exists at all, we find from (^9|) and ( p2| ) 

f = -1.3 x KT^Vf g A m-iM-iMpfT^. (93) 
In this case it is convenient to express the required conditions in terms of the inequalities on the temperature. 

(i) T< 2,.7c^g-§m§M*M^ v 

(ii) T« Osiycjs* ml Mi M*,. 

(iii) T<0.32g~£m*M%. 

(iv) T»1.2x 10 3 3- 4 mMMp/. 

(v) T « 0.22 fl J C ;^*m*M*M*,. 
(vii) r»2.3c^sr£m§M&. 

In addition, consistency of the assumption that stays in the time-dependent minimum (|8^ ) requires ry > r Q , namely, 

T > 5.0c^g"«miM'. (94) 
The e-folding number in this regime is expressed by an integral over the temperature as 

N S St. H T = ^c^g-'miM^M^ (tJ% T t ^) , (95) 

where Ti and Tf are the upper and the lower bounds on the temperature satisfying all the above inequalities, and 
we have set <j> ~ M. Inserting ( |94]) to ( |95|) we obtain iV <C 7.2 x lO^^g^^c^g^imi Mi M pl 2 . Now the consistency 
between (i) and (iv) imposes an upper bound on mM 2 as mM 2 <C 5.3 x 10~ 7 c^ 2 g a M Pl . From these two inequalities 

— 1 A. 33 

we find N <C 6.5 x 10 g^cFgis, which is much smaller than unity. 

Thus we do not see any possibility to implement the idea of warm inflation even in this model. 

V. CONCLUSION 

In the present paper we have examined feasibility of the warm inflation scenario from various view points. 

First we discussed how the viscosity term arises in the equation of motion of a scalar field in a thermal bath following 
Hosoya and Sakagami |l(| . Indeed such a term as C v <j) appears because it takes finite time for number density of 
particles interacting with <j> to relax to its thermally equilibrium value while <p is in motion. 

The viscosity term thus obtained could be very large at first glance, because it is not suppressed by any small 
coupling constants. As is evident in its derivation, however, this term appears as a result of a small correction to a 
sub-leading thermal correction term in the equation of motion. Since we know that even the leading thermal correction 
term does not lead to inflation, such a sub-leading viscosity term is not expected to play an important role and to 
yield an inflationary regime of a new type. 

If one neglects this fundamental feature and solves the overdamped equation of motion including the term C v <f>, 
one may find solutions indicating a possible emergence of warm inflation. However, we have found that such solutions 
violate the adiabatic condition that the scalar field should not change significantly in the relaxation time of the 
particles interacting with it. This condition is necessary for the derivation of the viscosity term C v (f>. Thus, the 
equation of motion incorporating the viscosity term C v <f> fails in the regime where it could describe warm inflation. 

If, on the other hand, we attempt to realize warm inflation in a manner fully consistent with the field theoretic 
derivation of the equation of motion, we inevitably find that the number of e-folds of inflation is constrained to be much 
smaller than unity, mainly due to the difficulty to satisfy the high-temperature condition and the adiabatic condition 
simultaneously. Even the shifted field model, which has been proposed to relax the high-temperature condition 
turned out to be no exception. 

Thus, our results as a whole show that it is extremely difficult or perhaps even impossible to realize the idea of 
warm inflation. 
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VI. APPENDIX 



In this appendix we derive the viscosity coefficient arising from Yukawa interaction. First we calculate the fermion 
self energy E^,, which is expressed as 



f d 3 k 

Ev,(p,r 2 - Ti) = f J -p^T3 S ( T 2 ~ n,k)G(T 2 - n,p-k), 



(96) 



in terms of temperature Green function of if), S(t, k) and that of 0, G(r, q). We use the following spectral represen- 
tation of Green functions. 



S(r,k) -- 
a(u>, k) = i 
G(r,q) -- 
p(w,q) = i 



du> 
2tt~ 



a(u>, k)e- UT [(1 - n F (u)) 6{t) - n F (uj)9(-T)} 



lujo — k-y + uj-f — k^ + 



(ui + ie) 2 — k — m 2 (uj — ie) 2 — k — m 2 
^p(u, q)e-" T [- (1 + n B (u)) 6(r) - u b (uj)6(-t)} 



(uj + ie) 2 — q 2 — m , (lu — ie) 2 — q 2 — m 2 



(97) 
(98) 
(99) 
(100) 



Inserting them to (|9£|) and applying the Fourier transform we find the self energy in Matsubara representation 



E^(w h p) = [ dre^E^r.p) 
Jo 



-r 



d 3 k f dcui du>2 



(2?r) 3 J 2tt 2tt 



a(ui,k)p(u! 2 ,p - k 



1 + n B (ui 2 ) - n F {u>{) 
iu>i - (u>i + UJ 2 ) 



(101) 



with uii = (21 + \)ttT . The above expression can be analytically continued to the retarded self energy E^(po + ie,p), 
and the imaginary part of the self energy is given by the discontinuity 



ImE^(p) = ~r [E^(p + ie,p) - E (p o - ie,p)] 



(102) 



Explicitly, we find 

ImE^(p) = tt/ 2 



d 3 k f dui\ dui 
(2 
d 3 k 



k I ^l^ ff ^ i)fc ^ W2i p_ k j ^ + n B (uj 2 ) - n F (uJx)) S (po -u>i- lo 2 ) 

7T)° J Z7T Z7T 

dw\dui 2 (wi7o - k~f + m^) (1 + n B {ui2) - n F (u)x)) sgaui\8 (ui\ - k 2 - m^) 



(2tt) 3 

xsgnw 2 <$ {oj\ - (P - kf - ml) 



(103) 



Putting p = (po,0,0,p3), \k\ = k m , and \p\ = p m (— P3), and in the limit and are negligible compared with 
p m , we find 



TT f 2 T 2 TT f 2 T 2 

ImE^(p) ~ — — — (70 - 73) = ^7—2- t> = r A 



64 pr, 



64 p 2 



(104) 



for p m = po < T. 

Then the dressed retarded Green function reads 
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S r (po,p) 



(1 - iT) ((p + «e)7° - Pi) - m 4:T ' 



(105) 



in momentum representation, where m^T = + Re£^ is the finite-temperature effective mass. The dressed spectral 
function is therefore given by 



a(p)= l (S R (p)-S^(p)) 

(1 - iT) fa + m^ T 



(1 + iT) ft + m^ T 



_ (1 - iT) 2 [(po + i£ )2 _ p 2] _ m 2 t (1 + , r) 2 [(p() _ Z£ )2 _ p2] _ m 2^ 



(106) 



= 2 



(1 - iT) jp + m^ T 
2u p {\ - it) 2 



1 



(1 + if) $ + m^T 



2<(i + ir) 2 Vpo-^; po + ^; 



with u>p — p m + irn^rpT /p m . The final expression applies in the limit \p\ = p m 3> ra^T- 

In terms of this spectral function the real-time finite-temperature Green function reads in momentum representation 



S F (t,p) = -i J ^.e-^Hl -n F (p ))9(t)-n F (p )9(-t)]a(p) 



(107) 



(l-if)K-yy) +ffH>r (l + iDK-pr)^ ( _ } K( 

2 Wp (i-ir) 2 v p7 2 Wp (i + ir) 2 v p;; 

(i- i f)K- P7) + ^ r (1 + i f)(- t , ; -p 7 ) + m ^ _ 

2w p (i-ir) 2 1 py 2 W ;(i + ir) 2 v ty p " 

Therefore the reading term to the viscosity due to Yukawa coupling is given by 

d 3 p 



9(t) 
<?(-*)■ 



2/> / dt'(t'-t) 



(2 _ )2 lmS^(t-t\p)S^(t'-t,p) 



2/ 



d 3 p 



2r(t'-t) 4w v»r 



(27T) 2 

64 i /*°° J 



df (f - t)e~ 



f,PPn 




ir 



'(e/3pm + l)^ T 3 



= ^C(3)T0~ 11.2T0, 



(108) 



where 



r = "4ta_ *f 2 T 2 ml T 



64p3 



(109) 
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